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ON THE SOLAR ROTATION RATE IN THE UPPER CONVECTION ZONE 
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ABSTRACT 



We present results on variation in rotation rate in 
the upper convection zone using data from GONG 
and MDI/SOHO covering a period of more than four 
years. We find that the first few odd-order splitting 
coefficients vary systematically with the solar cycle. 
The rotation rate near the solar surface calculated 
from analytical methods agrees well with that of in- 
version techniques. The residual rotation rate in the 
outer layers seem to be correlated with the solar ac- 
tivity. 
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1. INTRODUCTION 



With the advent of helioseismology and specially 
over last few years remarkable progress has been 
made in characterising the differential rotation of 
the sun. Analysis of the first 4 month's Global Os- 
cillation Ne twork Group (GONG) data (Thompson 
|st al., 199q ) have qualitatively confirmed the results 
of earlier analyses that the surface differential rota- 
tion changes to rigid rotation below the outer 30% of 
the solar radius. In addition, data from GONG net- 
work and medium-^ program of Michelson Doppler 
Imager (MDI) on board Solar and Heliospheric Ob- 
servatroy (SOHO) have allowed more detailed study 
of the rotational dynamics. This is achieved due to 
the splitting of the sun's global oscillation frequen- 
cies by large scale flo ws. Details o f this can be found 
in a recent review by Howe (1998). 
I 1 

Woodard & Libbrecht (1993) studied the time varia- 
tion of equatorial rotation rate for different £ ranges 
and suggested a small variation from year to year. 
Inversions for the rotation rate in solar interior have 
confirmed that the temporal variations in the ro- 
tation rate penetrate to some what deeper layers 

flbomrc et al., 2000 ; 
the changing pattern 
of solar rotation rate in the interior also agrees with 
the torsonal oscillations obs erved at the solar surface 
(Howard fc LaBonte, 1980). These zonal variations 



( Howe, Komm fe Hill, 2000 
Antia fe Hasu, 2000| ) lurthc 



by Kosovichev fc Schou (1997 ) using MDI /-mode 
data. Toomrc ct al. (2000) ) have confirmed these 
zonal flows using both p- and /-mode splittings de- 
rived from MDI data. 

In this paper, we analyse the odd order splitting coef- 
ficients for detecting any periodic or systematic vari- 
ation with s olar cycle. Ext ending the analytical for- 
mulation of Morrow (1988 ) to include coefficients up 
to 9th order, we infer the solar ratation rate as a 
function of depth and latitude. These rotation pro- 
files are compared with those obtained from a 1.5d 
inversion technique. We also estimate the variation 
in the rotation rate from the solar minimum to max- 
imum period of the current cycle. 



2. THE DATA AND ANALYSIS TECHNIQUE 



It is known that the solar differential rotation and 
other symmetry breaking factors like magnetic field 
can lift the degeneracy of the solar acoustic modes 
and split the eigen frequencies. Individual mode 
splittings can be represented by polynomial expan- 
sion 



s=l 



eP a (m), 



(1) 



where v n j is the mean multiplet frequency, P s (m)'s 
are orthogonal polynomials of degree I and the ex- 
pansion coefficients a s ,n,i are known as the splitting 
coefficients. The odd-order coefficients measure the 
solar rotation while the even-order coefficients probe 
the symmetry about the equator. The nonzero val- 
ues of these even coefficients reflect t he pole-equator 



asymmet ries in the solar structure ( Dziembowski et 
al., 2000|). 



of the sun's differential rotation were first observed 



We us e data from b oth GONG ( [Hill et al., 1996| ) and 
MDI ( |5chou, 1999 ) to determine the rotation rate in 
the outer convection zone. Each of the 45 GONG 
data sets covers a period of 108 days and is centered 
on dates 36 days apart. The month 1 starts from May 
7, 1995 and month 46 ends on December 23, 1999. 
The MDI data consists of twenty 72 days time series 
starting from May 1, 1996 and ending on August 31, 
2000 and includes the breakdown period of mid 1998. 
In case of the GONG data sets, splitting coefficients 
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Table 1. Correlation statistics for odd order coeffi- 
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are fitted up to ag while in MDI data, the fits obtain 
coefficients up to ois- In the present work, we further 
restrict the frequency range to 1.5 < v < 3 mHz and 
i > 70. 



3. RESULTS AND DISCUSSIONS 



3.1. Variations in the Splitting Coefficients 



L" 



Howe, Komm, & Hill (1999) studied the variation in 
GONG odd-order a coefficients but could not find 
significant correlation with solar activity measures . 
Later the same authors ( Howe, Komm, fc Hill, 2000 ) 
reported small trends up to ais coefficients. We have 
carried out a detailed analysis of these coefficients 
and observe that first few odd-order coefficients (up 
to 07) from MDI show significant correlation with 
activity measures (see Table-1). However, GONG 
odd-order coefficients reveal weaker correlation with 
activity. The temporal variation in odd-order coeffi- 
cients is shown in Figure 1. We notice that there is 
no clear trend in higher order coefficients while lower 
order coefficients vary in a systematic way. Fur- 
ther, the ct3 coefficients from MDI reveal strong anti- 
correlation with activity while these from GONG do 
not show any systematic trend. 



3.2. Solar Rotation Rate 



The solar rotation rate using the helioseismic data 
are generally obtained through two different meth- 
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Figure 1. The variation of mean odd-order splitting 
coefficients with time. The MDI data is represented 
by solid line while the dashed line shows the GONG 
coefficients. The la error bars in MDI data are not 
clearly visible since these are very small. 



ods: in forward approach the frequency splittings arc 
computed for a chosen solar rotation model and then 
compared with the observed splittings. In the inverse 
method the measured frequency splittings are used 
directly to produce a single function for the angular 
ve locity. In this s tudy, we use the analytical method 
of Morrow (1988) by including higher order splitting 
coefficients up to ag where the appropriate combi- 
nation of odd order splitting coefficients reflects the 
depth variation of angular velocity at chosen lati- 
tude 9 (see Tripathy, Jain, & Bhatnagar, 2000 for 
details). 



It is conventional to express the observed rotation 
rate in terms of even powers of cosO: 



(2) 



where 9 is the latitude and m is restricted to 2. With 
the availability of higher order a-coeffcients, we ex- 
tend and retain terms up to m — 4. In the asymptotic 
limit, we can express surface constants A4 in terms 
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Figure 2. The variation of average solar rotation rate 
calculated from five term fits in even power of cos 9 
are shown as a function of latitude. The line style 



have the following meaning: solid line for Snodgras 
long dashed line for GONG, dash- dot- dot- dot 
line for MDI and short dashed for GONG using 1.5d 
RLS inversion technique. 



of splitting coefficients as: 

A a = a\ + a 3 + a 5 + a-; + a 9 , 
A x = - [5a 3 + 14a 5 + 27a 7 + 44a 9 ], 
A 2 = 21a 5 + 99a 7 + 286a 9 , 
429 

A 3 = -[-—-07 + 57209], 
5 



2431 



-ag. 



(3) 
(4) 
(5) 

(6) 
(7) 



Various anttors have calculated ihese surface co- 

3 ffi 
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Brown et al 
56.7 nHz and 



cmcienLs using 

only three terms in the expression (2 
(19§|) found A = 462.8 nHz, A x = - 
A 2 = - 75.9 nHz for r > 0.723 R Q using Call K 
intensity data taken from South pole. Based on 100- 
day _obfiervay^nsmade at BBSO for I between 10 to 
60, [Libbrccht (1989| ) found the best fit with A = 461 
nHz, Ax = - 60.5 nHz and A 2 = - 75.4 nHz. We 
have derived these constants from the GONG data 
for the 70 < I < 150 and obtained the best fit with 
A = 460.55 nHz, Ax = - 55.76 nHz and A 2 = - 154 
nHz, A 3 = 220.36 nHz and A 4 = - 164.43 nHz. Simi- 
larly, MDI data yields A = 458.43 nHz, Ax = - 40.58 
nHz and A 2 = - 182.81 nHz, A 3 = 223.893 nHz and 
A4 = — 136.09 nHz. The values obtained for Aq are 
in close agreement with the earlier values. 

In Figure 2, we plot the average surface rotation rate 
as a function of latitude using the derived coeffi- 
cients. In the same plot, we have also shown the 
results from 1.5d RLS inversion an d Doppler sur- 
face measurements ( Bnodgrass, 1984 ). It is clear that 
the rotation rate changes significantly in mid latitude 
while the change in rotation rate near pole and equa- 
tor is small. We find that the rotation rate derived 



Figure 3. The temporal evolution of residual rotation 
rate near the surface at two different latitudes using 
both GONG (diamonds) and MDI data (stars). 



from the analytical approach agrees well with other 
results. However, the inverted rotation profile de- 
parts from other rotation rates beyond the latitude 
of 70 degree, probably due to the resolution limita- 
tion in inversion techniques. 

Figure 3 shows the residual surface rotation rate 
at two different latitudes obtained after subtract- 
ing the average angular velocity. It is evident that 
at higher latitudes, the residual rotation rate, com- 
monly known as zonal flows, is time dependent. This 
has a magnitude of approximately 3 nHz at equator. 
The residual rotation rate at equator for two di fferent 
depths is calculated u sing 1.5d RLS inversion ( Antia, 



Basu, fc Chitre, 200C) and is illustrated in Figure 4. 



It is clear that in the outer layers, the residual flows 
has a systematic variation while at deeper layers, the 
variation docs not appear to be systematic. Recently 
Antia fc Basu (2000[ ), |Howe, Komm fc Hill (2000| ) 
using GONG data, and tToomrc e t al. (2TOol] " using 



MDI data, have also reported a small but significant 
time variation in the rotation rate. 



4. CONCLUSION 



The first few odd-order coefficients are found to vary 
systematically with the solar cycle. We also detect 
a small but significant variation in the rotation rate 
derived from the linear combination of odd order co- 
efficients over a period of four and half years. The 
residual rotation rate in the outer layers seem to be 
correlated with the solar activity while no such sig- 
nature is found at the base of the convection zone. 
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Figure 4- The temporal evolution of the residual ro- 
tation rate at equator for two different radii as in- 
ferred from GONG data using 1.5d RLS inversion 
technique. 
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